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Abstract: Recent Holographic studies have shown that N=4 super Yang-Mills theory cou- 
pled to fundamental matter with finite chemical potential undergoes a first order phase tran- 
sition. In this paper, we study Nf D6 probe branes with or without electric field on it in the 
black D4 brane background compactified on a circle with supersymmetry breaking bound- 
ary condition. At energy scales much lower than the compactification scale, the dual gauge 
theory is effectively four dimensional non-supersymmetric SU(iV c ) Yang-Mills theory coupled 
to fundamental matter with or without baryon number charge. Within the supergravity ap- 
proximation, the decoupling of the Kaluza-Klein modes is not fully realized. For chemical 
potential fib < N c M q there is a line of first order phase transitions from stable meson phase 
to unstable meson phase. On the other hand for > N c M q there is no phase transition and 
mesons are unstable. A peculiar and interesting property of this system is that for a certain 
range of chemical potential fif, < N c M q , a new phase transition appears in the unstable meson 
phase. This phase transition is characterized by a discontinuous change of unstable meson 
lifetime. 
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1. Introduction 

Thermal phase structure of strongly coupled SU(iV c ) super Yang-Mills theory at finite chem- 
ical potential and finite temperature may be studied by using holographic duality [3, 4]. The 
holographic method provides a powerful framework to study a broad class of large N c , strongly 
coupled gauge theories with a small number Nf of fields in the fundamental representations 
[9, 10, 11, 12, 13, 14]. The gravity dual of these gauge fields appear as Nf probe Dq-branes 
on the near horizon geometries of iV c black Dp-branes. 

Particularly interesting physics appears when the system undergoes a confinement /de- 
confinement phase transition. Above this phase transition temperature Td ec , the gluons and 
the adjoint matters are deconfined and the dual geometry contains a black hole [5]. Since our 
main understanding of QCD largely comes from the behaviour of fundamental matter, it is 
interesting to ask how they behave around the deconfinement phase transition. The thermal 
properties of fundamental matter in the deconfinement phase can be extracted by studying 
probe D-branes in this black hole background. 

In [11, 12], it was show that at zero baryon density, fundamental matter undergoes a first 
order phase transition at Tf un . In the lower temperature phase, the probe D-brane sits entirely 
outside the black hole horizon, which is called a Minkowski embedding (see Figure 1), and the 
mass spectrum of meson in this phase is discrete and has a mass gap [11, 12]. In the higher 
temperature phase, a part of the probe D-brane fall through the black hole horizon, which is 
called a black hole embedding, and the mass spectrum of meson in this phase is continuous 
and gapless, characterized by quasinormal modes [12, 15, 16]. We emphasize that this is 
in the deconfinement phase. That means that after the gluons and the adjoint matters are 
deconfined, the fundamental matter is still in bound states at temperature T^ ec < T < Tf un . 
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The meson bound states in the deconfinement phase are also illustrated in related models 
[18]. On the other hand if Tt un is smaller than Td ec , the confinement /deconfinement phase 
transition and the meson melting take place simultaneously. This phenomenon, the stable 
meson bound states in a deconfinement phase Td ec < T < Tt un , is also found in lattice 
QCD[20], which suggests that bound states of heavy quarks survive after the deconfinement 
phase transition up to a few Tf un . 

Motivated by these agreements with 
QCD, other region of the QCD phase di- 
agram, i.e., finite temperature with fi- 
nite baryon density nj, or finite chemical 
potential fi q = fJ,b/N c was investigated 
in [2, 21, 22, 23, 24, 25, 26], see also 
[19]. The introduction of a finite baryon 
density corresponds to the presence of 
a U(l) gauge field, which is a diagonal 
part of the gauge group U(N/), on the 
probe D-brane. 

In [2], this direction is explicitly in- 
vestigated in the D3/D7 system. The 
boundary gauge theory is SU(A r c ) A/"=4 
super Yang-Mills theory with Nf J\f=2 
hypermultiplets in the fundamental rep- 
resentation. There, physical quantities 
such as quark condensate in finite n& con- 
tinuously change from those in n& = 0. 
One of the most crucial difference be- 
tween zero and nonzero baryon density 
systems is in the presence finite baryon density, Minkowski embeddings are unphysical and 
do not play any role, while black hole embeddings cover the whole range of the temperature 
above deconfinement phase transition. This is different from n& = case where black hole 
embeddings cover only high temperature region and Minkowski embeddings cover only low 
temperature region [11, 12]. In the overlapping region, there is a phase transitions from a 
Minkowski embedding to a black hole embedding. 

The reason why Minkowski embeddings are unphysical in the presence of finite baryon 
density is as follows; The electric field on the D-brane represents dissolved fundamental 
strings. Since Wess-Zumino coupling is inactive in this brane configuration, fundamental 
strings never 'leak' from the brane and the local baryon number density is still nonzero above 
the horizon. Since the fundamental strings cannot terminate, we have to have fundamental 
strings stretching from the probe D-brane to the horizon. At this junction point, the tension 
of the probe D-brane is always smaller than that of the fundamental strings. That means that 
there is no force balanced configurations in the Minkowski embeddings at finite baryon den- 
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Figure 1: Probe D-brane configurations in the pres- 
ence of the black hole. Gray circle represents the black 
hole horizon. For the black hole embeddings, tempera- 
ture decreases from the bottom to top. The dark green 
colour represents Minkowski embedding which is en- 
tirely outside the horizon. At lower temperature, the 
probe D-brane is entirely outside the black hole hori- 
zon, i.e., the Minkowski embedding. (dark green line). 
At higher temperature, the probe D-brane partially fall 
through the black hole horizon, i.e., the black hole em- 
bedding. 
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sity and the brane must fall down through the horizon. See [2, 1] for more detail discussions. 
Instead at very low temperature a black hole embedding mimics a Minkowski embedding, i.e., 
a very long and narrow spike, which corresponds to a bundle of dissolved fundamental strings 
on the probe D-brane, stretches down to the horizon. For smaller rib, there is a first order 
phase transition from a black hole to another black hole embedding, which is very similar to 
a Minkowski to a black hole embedding in the zero baryon density case. Above the critical 
density n£, there is no phase transition anymore. 

However it was shown in [2] that there is an unstable region near the phase transition 
line. It is expected that this unstable region finally decays to other stable state. However 
since this region minimizes the free energy, there is no state to decay and the final state was 
missing. 

This puzzle was recently addressed in [1]. In the grand canonical ensemble, we fix not 
a baryon density but a chemical potential. Since a Minkowski embedding can have a finite 
chemical potential without a baryon density, this embedding is physical in the grand canonical 
ensemble and plays an important role. The black hole embeddings cover only larger chemical 
potential or higher temperature region. On the other hand, the Minkowski embeddings cover 
the whole value of fi with temperature lower than Tf un (see fig. 2 in [1]). There are phase 
transitions inside the overlapping region from Minkowski to black hole embeddings. One of 
the most remarkable results is that the unstable region is thermodynamically unfavourable 
in the grand canonical ensemble. Since in the thermodynamical limit, the grand canonical 
and the canonical ensemble should give the same answer [28], the unstable region is not a true 
ground state but it should be replaced by an inhomogeneous phase of a stable black hole 
embedding and a Minkowski embedding in the canonical ensemble. 

The aim of this paper is to investigate the generality of the phase structure in Dp/Dq 
system. We analyze the D4/D6 system. Physical properties such as the phase transitions from 
a black hole to a black hole embedding in the canonical ensemble and from a Minkowski to a 
black hole embedding in the grand canonical ensemble and the existence of the unstable region 
are the same. However the inhomogeneous phase has a more complicated and interesting 
structure. There is a new additional phase transition. Figure 2 shows the phase diagram of 
the D4/D6 system. As in the case of the D3/D7 system, the black hole embeddings cover 
only high temperature or high chemical potential region. The boundary of the black hole 
embeddings is plotted in the green line. In the region surrounded by the blue line, there are 
more than one value of n q for each ji q /M q and T/M. The red line is the line of the phase 
transitions. On the scale of (a), the difference between the red and the green line is very subtle. 
An interesting structure is found in (b) and (c). Near n q /M q = 0.16, T/M = 0.77, the red line 
separates into two branches. The upper branch represents phase transitions from black hole 
to black hole embeddings. In the field theory side, this corresponds to a discontinuous change 
of meson lifetime. Below this phase transition temperature, the spectral function of mesons 
would show sharp peaks representing longer lifetime quasiparticles, while above the phase 
transition, the peaks would be flatter and the lifetime of the quasiparticles would be shorter. 
The lower branch represents phase transitions from Minkowski to black hole embeddings. We 
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will discuss the relation between these phases and the unstable region in Section 3. 
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Figure 2: Phase diagram: Quark chemical potential n q /M q versus temperature T/M. The red line 
separates the phase of Minkowski embeddings (small temperatures, small fi q /M q ) from black hole 
embeddings. Figure (b) zooms in on the region near the end of this line and also depicts the boundary 
of the region accessed by the black hole embeddings (green) and a small region (enclosed by the blue 
curve) where more than one black hole embedding in available for a given value of fi q and T. Figure 
(c) zooms in on the region near the branches of the red lines. The two branches of the red lines show 
the phase transitions, the lower one is from the Minkowski embeddings to the black hole embeddings 
and the upper one is from the black hole embeddings to the black hole embeddings. 



2. Holographic framework 

In this section we explain the D4/D6 system. Closed string modes on the near horizon 
geometry of D4 brane is dual to maximally supersymmetric five dimensional super Yang- 
Mills theory with large N c and at strong coupling on the boundary . In order to obtain four 
dimensional gauge theory at low energy scale, we compactify one space direction, say X4, with 
a radius M^} R . For energy scale E « Mkk, the boundary gauge theory is effectively four 
dimensional. Since we want to study the thermal properties of fundamental fields coupling to 
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Yang-Mills field, we choose the supersymmetry-breaking spin structure, i.e., the anti-periodic 
boundary condition for fermions in X4 direction. The fermions in the vector multiplet get 
masses of order Mkk at tree level and the scalars possibly get masses of order g\N c MKK at 
one loop level. So the gauge theory is effectively non-supersymmetric four dimensional Yang- 
Mills theory at low energy scale. 1 By fixing the boundary geometry, this system undergoes a 
phase transition from a thermal AdS phase to an AdS black hole phase at certain temperature 
Tdec- In the dual language, this transition is deconfinement /confinement phase transition [5]. 
This system describes only adjoint matters. To approach QCD in this holographic framework, 
quarks or fundamental matter must be included. The introduction of fundamental matter 
in the holographic framework is demonstrated in [8]. According to their work, a probe D- 
brane plays a role of fundamental matter in the dual picture. In our case, inserting Nj D6 
branes into this black-D4 background corresponds to the coupling of Nj hypermultiplets in 
the fundamental representation to the Yang-Mills fields with the gauge group SXJ(N C ). The 
quarks arise from the lightest mode of the fundamental strings connecting between D4 and 
D6 branes. In the decoupling limit, the mesons are dual to the fundamental strings whose 
both ends are attaching on D6 branes. This means that the fluctuations on the D6 branes 
describe the mesons in the dual gauge theory. 

2.1 Black D4 background 

The supergravity solution corresponding to the decoupling limit of N c coincident D4-branes 
in the black hole phase is in the string frame 



ds 2 



3/2 



f 2 

■-jdt 2 + fdx\ 



(!) 



+ 



r\3/2 ?i/ 3 

7) f ^W + ^ 



where / = 1 - |f and / = 1 + |± 2 . 

uq is the location of the horizon and the AdS radius L is given by 

L 3 = g s N c TTl 3 s . 



(2.1) 
(2.2) 



(2.3) 



Hawking radiation appears in this background with the temperature given by the surface 
gravity T = n/2ir, or by the regularity of the Euclidean section, 



T = 



4vr V L 3 ' 



(2.4) 



1 However the Kaluza-Klein modes do not decouple within the supergravity approximation [5]. For example, 
the lightest glueball spectrum is of the same order as the strong coupling scale [7]. 

2 This metric (isotropic coordinate) is related to the ordinary coordinate as («op) 3 ^ 2 — u z ^ 2 = \J « 3 — Uq 
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This temperature is identified with the temperature of the boundary gauge theory. The 
dual boundary gauge theory is maximally supersymmetric Yang-Mills theory on R 1,4 with 
coordinates {t, X4}. 

As we explained above, we compactify on a circle with radius 1/Mkk with supersymmetry 
breaking boundary condition. At low energy E « Mkk the effective field theory is four 
dimensional non-supersymmetric Yang-Mills theory. The string coupling constant in the bulk 
theory and the gauge coupling constant in the dual boundary theory is related through 

gl = (2vr) Vs, A = &M KK (2.5) 

where g s = e^ 00 and #4,5 are the Yang-Mills coupling constants in four and five dimensions. 

2.2 D6-brane embeddings 

We introduce Nf D6-branes in this background in the following configuration; 
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We split the AdS radial direction and S* 4 into two parts so that SO (3) x SO (2) symmetry 
is manifest. We define the radial coordinate r in (x 5 ,x 6 ,x 7 ) directions and R in (x 8 ,x 9 ) 
directions where p = q/uq is a dimensionless coordinate. 

p 2 = r 2 + R 2 , r = psm6, R = pcos6, (2.6) 



and 



dp z + p z dfli = dp' + p\dQ A + sin 2 6dW 2 + cos 2 6d^ 



dr 2 + r 2 dQ^ + dR 2 + R 2 dd> 2 



(2.7) 
(2.8) 



The probe D6 brane embedding is parametrized by a function, x{o) = cos #(£>). The induced 
D6 metric is 



ds' 



2 1 (qV 2 



2 \L 



-Ldt 2 + fdx 2 



,X 3/2/1/3 
' [ Q J 2V3 



Q 2 X 2 



i + ~ ~~7 ) dg 2 + g 2 (l- x 2 ) dn 2 

(2.9) 

Now we introduce U(l) gauge field A(g)dt on the probe D6 branes. The D6 brane action 
per unit spacetime volume of the gauge theory is 



l m = -N f T 6 J e-*^/-det(P[G} ab + 2Tra>F ab ) 



-NfTe^ J dgf^g 2 (l - X 2 ) J f*f-V* (l + - 



Q 2 i&_ 

r 



-kA 2 



(2.10) 
(2.11) 
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where 

k = (2 2 / 3 2tt/ 2 ) 2 , (2.12) 

and the D6 brane tension is 

T ' = A- <2 ' 13) 

Let us analyze the behaviour of the electric field A. The electric displacement d is defined 

by 

d N f T 6 Q 2 -X 2 )kA (214) 

4 ^/2/-2/3 (1 + £g) - kA> 

The electric displacement represents the number density of dissolved fundamental strings 
in the D6 branes. With the similar analysis to [2, 1] the electric displacement identified with 
the string number density n q : d = n q . Since the action does not explicitly depend on A, the 
equation of motion for A is d =constant, or 



A = V x (2.15) 

^fc+(2^)Vv(i-x 2 ) 2 

(2.16) 

Asymptotically x g° es to zero and / and / go to one. The leading order of the equation is 
then 

A = N . T L k - + o(^] (2.17) 



N f T 6 n 2 k 2 ' \ p 3 

4 * 



So the asymptotic behaviour of the gauge field is 



where \i is 



4d 1 . 
A~v q - — —— - + ■■■ 2.18 

NfT 6 U 2 k g 



d / dn V 1 X (219) 

Juo l d2k + j^A^WX 2 /8/3 ^ 4 (1 _ x2) 2 



Here we set A(uq) = 0. This condition comes from the regularity of the one form at the 
horizon. The horizon contains a bifurcation surface where a Killing vector dt vanishes. In 
order for the gauge field, which is a one form A t dt, to be well defined, the component A t must 
vanish there. See also [2, 12]. 

According to the holographic dictionary, the dual operator is schematically 

O q = ^ty + q^tq (2.20) 
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where ip and q are left and right Weyl fermions and scalar fields (see also [2]). From (2.18) 
and (2.20), we identify p q with a quark chemical potential. 

For convenience, we also use dimensionless quantities and a gauge field, 

7 4 d \fk ~ s/k . 

d = fi = —ii, A = — A. 2.21 

N f T 6 n 2 Vkul u n 



The equation of motion for x is 




(2.22) 



kA 2 \ 

- 2 K - = (2.23) 

f 2 f -2/3 ) 

The boundary conditions for x{p) a t the horizon is determined by the regularity and it gives 
xlp—i = Xo an( i dx/dp\ f) _ l = for < x < 1- For Minkowski embeddings it is convenient 
to use R, r coordinates instead of %, p coordinates. The boundary conditions for Minkowski 
embeddings at r = are -R| r=0 = Ro an d dR/dr\ r=0 = for 1 < Rq. The asymptotic forms 
(p — > oo) for x is 

X-^ + j^ + O^) (2.24) 

= f + # + °(>)> ( 2 - 25 ) 

where we define dimensionless quantities 

m 3 2 m _ c 3 4 c 

m = ^ = (4vr) 2 L3T 2 ' C = ^ = (4tt) 4 L6T 4 - ( } 

Holography relates these quantities to a quark mass and a condensate [10] by 

M, = J?2L (2^27) 

23 7T/ 2 

<O m > = -2 5 / 3 TT 2 l 2 s N f T m u 2 c (2.28) 

A bare quark mass M q is an asymptotic distance between D4 and D6 brane in flat space. 
The operator O m is the variation of the mass term in the microscopic Lagrangian, i.e., O m = 
—dM q £, and the schematic form is 

O m = i>4> + q^q + M q qU, (2.29) 

where $ is one of the adjoint scalars. See also [2, 12]. We assume that the vacuum expectation 
value of the fundamental scalar fields q vanishes because whenever the scalar fields acquire 
the nonzero expectation value, the energy density increases [10]. In that case, < O m > is equal 
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to < x/jip > and represents a quark condensation. Given the above result, we show various 
figures in terms of T/M = l/\ffti. 

Let us study more on the chemical potential. Once the equation of motion for x{o) 
(2.23) is solved, the chemical potential is obtained from (2.19). In general, we have to resort 
to numerical calculation. However, we can extract analytic properties in the limiting cases 
of small and high temperatures. First we consider low temperature or large bare quark mass 
limit T/M — > 0. As shown in Figure 1, at very low temperature the probe D-brane goes up 
straight from the horizon uq to ~ m. The main contribution to the distance between the D4 
and D6 brane, or T/M, comes from this part. Under this approximation, 



/•oo 

d / dg 

J UQ 



//-Ia/I + jlx. 



d?k + 



~ d / dg- 

J u 



m 



d?k 
/Vk = M q 



N f T 6 Q 2 k 
4 



/ 8 /V(i-x 2 ) 2 



(2.30) 

(2.31) 
(2.32) 



Regardless of d, fig goes to M q . This is consistent with Figure 3 where all curves with different 
d meet on the vertical axis at \x q /M q = 1. On the opposite limit, i.e., high temperature or 




Figure 3: Chemical potential n q /M q versus T/M for various values of d, increasing from bottom up: 
d= 10- 2 ,1,10, 100, 10 4 . 



small bare quark mass limit T/M — > oo, the probe D6 branes are almost flat and intersect 
the horizon at the equator. Plugging x — an d uq >> 1 into the equation, the leading 
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contribution is 



**v^*k + °w (2 - 33) 

- 9 nq +0{nl) (2.34) 



4N f N c T 2 

There are some notable points. The first point is that for any value of p q goes to zero 
as temperature goes to infinity. In this sense, Figure 3 might be misleading. Each line in 
the figure shows p q /M q versus T/M for fixed d not d. In each line, higher temperature 
corresponds to a higher baryon number density and lower temperature corresponds to a lower 
baryon number density. The second point is for fixed p q , rib changes as ~ T 2 . Since in the 
grand canonical ensemble we survey fix \i q planes, this equation suggests that the baryon 
density monotonically increases at very high temperature for any fixed ii q /M q . 

3. Thermodynamics 

We move on to the thermodynamics of the D6 brane, or the fundamental matter. As we 
will see the canonical ensemble is not a suitable one for studying the phase diagram since it 
includes an unstable region. Our main focus here is the grand canonical ensemble. A similar 
unstable region can also be found in the D3/D7 system [1]. We will briefly mention on the 
canonical ensemble. 

Euclidean path integral of D6 branes give thermal partition function [27]. Since the clas- 
sical solution of the equation of motion is the saddle point of the path integral, the on-shell 
action gives main contribution to the Gibbs free energy W, i.e., W = TIe- From the grand 
canonical point of view, a brane configuration which minimizes the Gibbs free energy is 
thermodynamically favourable. However, naive calculation of thermal quantities include di- 
vergences from IR region. For example, we can clearly see that the action diverges as 



/ 



dpp 2 ~ %21. (3.1) 



In order to obtain physically meaningful, finite quantities, we need to renormalize them. 
The renormalization procedure of probe D-brane is studied in [29]. According to that, renor- 
malization is done by introducing boundary terms for probe branes on the cut-off plane so 
that the divergences cancel. For the D4/D6 case, this boundary term is concretely studied in 
[12]. Before starting the detail analysis, we define a normalization constant 

N = ^N f T m ul (3.2) 
Inserting the asymptotic expansion of x an d A into (2.11), the action is 

W*4-")( 1 + f'''-3- (3.3, 
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The contribution to the action from the gauge field is the term proportional to d 2 . This gives 
no new divergence. So we can use the same boundary term as that of d=0 case to renormalize 
the divergence. We apply the following boundary term per unit spacetime volume of the gauge 
field [12], 



h 



bound 



3 

ulQ 2 T 6 N f 



12 



where 7 is the induced metric at p = p v 



P 3 

rmax 



-jfh p max - 3mc 



,2 1 fU p max \3/2 / f(p max ) 2 „ ~ „ 

ds i=0\ T J I "77 \ +f{Pmax)dx 3 

\ J (Pmax) 



The total action is then, 
hot 

M 



heg hound 



G{m,A)-±(p* 



™2 . 



3mc 



(3.4) 
(3.5) 

(3.6) 

(3.7) 
(3.8) 



where G(rh, A) is 



G(m, A) I dp W (1 - x 2 ) y Z 2 /- 2 / 3 (l + f±p ) - l 2 - ^ 2 + 



2V/2 



m 
T 



(3.9) 



From the thermodynamical point of view, this action is identified the Gibbs free energy 
W(T, /j, q ) via W = TI E . 

In the canonical ensemble with fixed n q , we use the Helmholtz free energy. Similar to 
[2, 30], the Helmholtz free energy is associated with the Legendre transform of Ie- 



Ie_ _ Idq J dA 



(3.10) 



which is function of the temperature and the baryon density. We identify F(T, n q ) = TIe 
where F(T, n q ) is the Helmholtz free energy. Since there is no contribution from the electric 
displacement to the boundary term , the Legendre transformed action is 



Ie 
M 



= H(m,d) 



If 3 

rmin 



m 2 p m in ~ 3mc 



(3.11) 



where H(m, d) is 



/oo 
min 



dp 



ff 



'1/3, 



'1 + 



P 2 X 2 



~ 2 

m 



2V /8/3 /9 4 (1 _ x 2 ) 2 + j2_ p 2 + 



(3.12) 
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We evaluated the free energy numerically. The qualitative feature is very similar to that of 
the D3/D7 system (see Figure 9 and 10 in [2]). For smaller value of d/m 2 , there is a line 
of first order phase transition from black hole to black hole embeddings and above a critical 
value ^d/m 2 ^j there is no phase transition. Note that since d depends on a temperature 

n q /T A , we use a temperature independent quantity d/fh? ~ n q /M 2 as a parameter. The 
blue line in Figure 9 shows the phase transition line in the canonical ensemble. The critical 
temperature and baryon density is T/M = 0.7786 and d/fn 2 = 0.041. 

3.1 Phase transitions in the grand canonical ensemble 

As we mentioned in Section 1, the black hole embeddings cover only above the green line in 
Figure 2 while the Minkowski embeddings cover whole value of n q /M q below ~ Tf un . This 
suggests that there is a line of first order phase transitions for fj, q /M q < 1 in the grand 
canonical ensemble. Thermodynamically favourable embeddings minimize the Gibbs free 
energy for give temperature and chemical potential. In this section, we calculate the free 
energy for various n q /M q . 
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Figure 4: Free energy versus temperature for fj, q /M q = 0.14. The blue dotted (red solid) line 
represents the Minkowski (black hole) branch. The vertical line marks the temperature of the phase 
transition. 

Figure 4(a) shows the Gibbs free energy W normalized by M versus temperature T/M 
for jj, q = 0.14 in a broad view and Figure 4(b) shows a zoomed in around the phase transition. 
Here the normalization is N = MM. The red line is the black hole embedding and the dashed 
blue line is the Minkowski embedding. The starting point of the red line (T/M = 0.776, 
W/N = —0.0077) is d = 10~ 3 . As d becomes larger, the temperature becomes higher and the 
free energy becomes smaller, then at d = 5 X 10 -2 the temperature starts to decrease and the 
free energy starts increase. And at d = 0.17, the temperature turns to increase and the free 
energy turns decrease again. As fi increases, the starting point of the free energy of the black 
hole embedding at very small d goes in the left above direction and the three fold structure 
starts to form a swallow tail shape. At a critical value of /i g /Af g =0.15, the crossing point of 
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Figure 5: Free energy versus temperature for (a)/i g /M g =0.165 and (b)/i g /M g =0.175. (a) shows that 
there are two phase transitions. One is from the Minkowski embedding to the black hole embedding 
at T/M = 0.765 and the other is from the black hole embedding to the black hole embedding at 
T/M = 0.767. (b) shows that there is only one phase transition from the Minkowski embedding to 
the black hole embedding at T/M — 0.76. The vertical lines mark the temperature of the phase 
transitions. 



the swallow tail goes down below the line of the Minkowski free energy. Above this critical 
value, there are two phase transitions. The first one is from a Minkowski embedding to a 
black hole embedding. The second one is from a black hole embedding to another black hole 
embedding. A representative of this phase is shown in Figure 5(a). At [i q jM q = 0.165, there 
is phase transition from a Minkowski embedding to a black hole embedding at T/M = 0.765 
and there is a phase transition from a black hole embedding to a black hole embedding at 
T/M = 0.767. As [i becomes larger, the swallow tail shrinks smaller and smaller and finally 
at n q /M q = 0.175 the black hole to black hole phase transition disappears. The free energy 
at this value of fJ- q /M q is shown in Figure 5(b). There is a phase transition at T/M = 0.76 
and this is the only phase transition in this phase. These three fold structures in the D4/D6 
system should be compared to the Figure 6 which shows the typical phase transition for 
[i q < M q in the D3/D7 system (see also Figure 4 in [1]). The only phase transition is from 
a Minkowski embedding to a black hole embedding and there is no phase transition from a 
black hole to a black hole embedding. 

This feature can also be seen from other perspective. Figure 7(a) shows d/fh? versus T/M 
diagram near the phase transition for fi q /M q = 0.16. As explained above, below a certain 
temperature there are only Minkowski embeddings. Around T/M = 0.7668 black hole em- 
beddings appear. However the Minkowski embeddings are still favoured. At T/M = 0.7673, 
the free energy of the black hole embedding becomes equal to the Minkowski embedding and 
the system jumps from d/fh? = to d/fh? = 0.0034 (A) As temperature increases, the system 
goes along the red line until it meets another phase transition point(B). At T/M = 0.7685, 
d/fh? = 0.0146 (B), the system jumps to a configuration which has the same chemical poten- 
tial but larger value of d/fh? = 0.0708(D) in a black hole embedding. Above this temperature, 
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there is no phase transition and as T/M becomes larger d/rh 2 becomes larger. 

We also show a d/m 3 versus T/M 
diagram for fi q /M q = 0.005 in the D3/D7 
system in Figure 7(b). As above, be- -o.did4 
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the free energy of the Minkowski embed- 
ding and the black hole embedding co- 
incides and the first order phase transi- 
tion occurs from d/rh s = to d/rh 3 = 
0.00073(G). Clearly [i q /M q is a two- valued 
function of T/M while it is a triple- valued Figure 6: Free energy versus temperature for fi q /M g = 
function in the case of the D4/D6 sys- 0.005 in the D3/D7 system. Fhere is only one turnover 
tern, and the region corresponding to (A)- in the black hole embeddings. Fhis should be compared 
(B) in Figure 7(a) is missing. with the D4/D6 system (Figure 4) where there are two 
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Figure 7: Baryon density versus temperature for (a) fi q /M q = 0.16 in the D4/D6 system and (b) 
fiq/Mg — 0.005 in the D3/D7 system. In the D4/D6 system, baryon density is a triple- valued function 
of temperature while in the D3/D7 system, baryon density is a double- valued function of temperature. 

In the rest of this section, we address another interesting property, i.e., the instability 
in the canonical ensemble, and combine it with the phase transitions discussed above. The 
conditions for the stability of the system are given 



dS_ 
df 



>0, 



dn b 
dn q 



T 



> 0. 



(3.13) 



The unstable configuration can be found by examining fi q /M q versus T/M and d/rh 2 dia- 
gram(Figure 8). This figure covers only a small range of the surface near the phase transition 
in the canonical ensemble. For fixed n&, or d/rh 2 , the system jumps from the top to the 
bottom of the surface at certain temperature. Thermo dynamical point of view, this means 
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Figure 8: Three-dimensional plot of the chemical potential, the temperature and the charge density 
determined by black hole embeddings. 

that only the region of the top of this fold with temperature lower the phase transition and of 
the bottom of this fold with temperature higher than the phase transition is favoured. Other 
region are thermodynamically unfavoured do not come into the stability discussion. An in- 
teresting point is that on the top of the fold and just below the phase transition temperature 
we can actually find a region where thermodynamically favourable but electrically unstable. 
This region should play a central role in the canonical ensemble. In [1], the similar problem 
was addressed. The resolution of the problem in the D3/D7 system was that minimizing 
the free energy in the grand canonical ensemble always picks out either a stable black hole 
embedding or Minkowski embedding. Hence the system never suffer from the instability. The 
question is whether the same is true for the D4/D6 system. After explaining Figure 9, we 
come back to this question. 

Figure 9 shows the black hole to Minkowski phase transition, dfiq/drib = 0, and the 
canonical phase transition line. The blue line is the canonical phase transition line. For 
fixed value of d, there is a black hole to black hole phase transition at this line. The green 
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Figure 9: Baryon density at the phase transition in the grand canonical ensemble (red) The blue line 
of phase transitions identified in the canonical ensemble at fixed n\,. The region enclosed by the green 
and blue curves corresponds to the unstable region in which (d/i q /dn q )T < 0. 

line is a part of the = line. The region surrounded by green and blue line is the 

thermodynamically favourable but electromagnetically unstable region. The red line is the 
mapping of the phase transition line in Figure 2 to d/fh? — T/M plane. Let us explain this 
line more detail. In Figure 4(b), there is one phase transition from the black hole embedding 
to the Minkowski embedding at T = 0.7742 and d/m 2 = 0.069. This corresponds to the 
red line between A and B in Figure 9. In this region there is only one value of d/m 2 for 
fixed T/M and at this point the phase transition occurs from a Minkowski embedding to a 
black hole embedding, i.e., from a point at d/m 2 =0 to a point on the red line with the same 
temperature. 

In Figure 5(a), there are two phase transitions. One is from the Minkowski embedding 
to the black hole embedding at T/M = 0.7649 and d/m 2 = 0.0034. The other one is from 
the black hole embedding at T/M = 0.76703 and d/m 2 = 0.0285 to the other black hole 
embedding at T/M = 0.76703 and d/m 2 = 0.07014 This configuration corresponds to a 
region between the two vertical black lines in Figure 9. The first phase transition is from a 
Minkowski embedding d/m 2 =0 to a point on the red line between E and G. Then the second 
phase transition occurs at a different temperature from a point on the red line between E and 
F to a point on the red line between F and B. 

In Figure 5(b), there is only one phase transition from the Minkowski embedding to the 
black hole embedding at T/M = 0.76 and d/m 2 =0.003. This corresponds to the left part of 
G on the red line. 

As in the D3/D7 system [1], the whole unstable region is surrounded by the red line. Hence 
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the grand canonical ensemble picks out either stable black hole embedding or Minkowski 
embedding; the thermodynamically favourable phases discussed above are all stable. The 
unstable region appeared in the canonical ensemble is misidentified the ground state since 
we restricted our analysis to the homogeneous configurations. The true ground state is an 
inhomogeneous phase of stable phases [1]. 

The true ground state of the region surrounded by the red curve A-B and the d/rh 2 = 
axis is an inhomogeneous phase of stable black hole (on the red curve A-B) embedding and 
Minkowski embedding (on the d/rh 2 = axis). The true ground state of the region surrounded 
by B-F-E-B is an inhomogeneous phase of stable black hole embedding (on the red curve B-F) 
and another stable black hole embedding (on the red curve F-E). The true ground state of 
the region surrounded by the red curve E-G and the d/rh 2 = axis, is an inhomogeneous 
phase of stable black hole embedding (on the red curve E-F) and a Minkowski phase (on the 
d/fh 2 = axis). Hence both in the grand canonical ensemble and the canonical ensemble, 
the unstable region is thermodynamically unfavourable. 

4. Conclusion and Discussion 

We have seen that there are several common physical properties in the Dp/Dq systems with 
finite charge density and chemical potential. In the grand canonical ensemble there is a 
line of the first order phase transitions from Minkowski to black hole embeddings found in 
[24, 25, 26, 1]. The black hole embeddings cover high T/M or high ii q /M q region while the 
Minkowski embeddings cover the whole ti q /M q region with T/M smaller than ~ Tf un . In 
the supersymmetric limit T/M — ► 0, fi q /M q always goes to one for any value of n&, i.e., the 
energy which is necessary to add one quark in the system is equal to its bare mass. In this 
limit black hole embeddings exactly look like Minkowski embeddings. There is a very thin 
and long spike stretching from Dq branes down to Dp branes. In the dual field theory point 
of view [12, 15, 16], the Minkowski embeddings correspond to the stable meson phase and he 
spectral function consists of a series of delta- function-like peaks, i.e., resonances centred on 
mass eigenvalues. On the other hand, the black hole embeddings correspond to the unstable 
meson phase and the spectral function is continuous. In the canonical ensemble, there is 
a line of first order phase transitions from black hole to black hole embeddings for charge 
density less than n£ and above this critical density there is no phase transition. Just below 
the phase transition temperature, there is an electrodynamically unstable region. However 
this unstable region is not the true ground state and should be replaced by an inhomogeneous 
phase. 

We also have seen that there is a difference between the D3/D7 and D4/D6 system. In 
the D3/D7 system for fixed \i q < M q , there is only one phase transition, from a Minkowski 
to a black hole embedding. On the other hand, in the D4/D6 system for certain range of 
fi q < M q , there are two phase transitions, one is from a Minkowski to a black hole embedding 
and the other one is from a black hole to another black hole embedding. Following the inves- 
tigation of [15, 17, 23], we can give a physical interpretation of this black hole to black hole 
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embedding phase transition in the dual gauge theory. The feature of the spectral function is 
characterized by the poles of the corresponding retarded correlators in the complex frequency 
plane. In the Minkowski embeddings, the poles are on the real axis and the spectral function 
is, as mentioned above, a series of delta-function-like peaks. In the black hole embeddings, 
the poles are located apart from the real axis. At lower temperature, the imaginary part of 
the poles are close to the real axis and the spectral function exhibits distinct peaks. As tem- 
perature increases, the poles move away from the real axis and the spectral function becomes 
featureless. Hence our black hole to black hole phase transition from lower temperature to 
higher temperature would correspond to poles jumping from closer locations to the real axis 
to farther locations. The peaks of the spectral function would become lower and the life time 
of the quasiparticles would become shorter. It would be very interesting to investigate this 
point concretely. Lattice QCD studies [20] suggest that meson bound states survive the de- 
confinement phase transition and the bound states dissolve at Tf un . As explained, Minkowski 
to black hole phase transition represents this meson dissolving phase transition. A very in- 
teresting question is then what the phase transition in QCD corresponding to the black hole 
to black hole phase transition is. It would be an exotic phase transition if it exists in QCD. 

Before ending the section, we comment one more common property in the D3/D7 and 
the D4/D6 systems. As mentioned above, in the T/M — ► limit, the chemical potential \x q 
goes to M q . This result is reasonable in the case of free fermions with zero density. However 
since our system is strongly interacting system with bosons and fermions, the reason is not 
very clear. In fact, the chemical potential is very close to the constituent mass at the phase 
transition points not only at T/M = but also for a broad range of temperature except the 
canonical phase transition region. Hence the quarks behave as if they were 'free' particles. We 
have put 'free' om quotes because the constituent mass takes into account all the corrections 
from the quarks and the adjoint plasma. The constituent quark mass is defined by the energy 
of a fundamental string stretching from a probe D-brane to the horizon in a Minkowski 
embedding. This constituent quark mass at T/M = on Minkowski embeddings with n q = 
was calculated in [12, 31, 1] The induced metric of the fundamental string stretching from a 
probe D-brane to the horizon in R direction is 

2 1 (u R\ 3/2 \ ^ ( L V' 2 /V3 

ds= 2{—) r7 dt \ + {^R~) ¥^ u ° dR - 

The Nambu-Goto action for the fundamental string is 

I fRo 

/ = -= / dRdt^-detP[G] 

27rZ| J 1 

-4/^(( 1 + 4) 2/3% - 22,s ) 

Identifying the constituent mass to minus the action per unit time of this static configuration, 
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we have 
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Figure 10: Comparison of the ratio M c /M q (blue) with the ratio yu q jM q (red) at which the phase 
transition from a Minkowski to a black hole embedding takes place. The two curves essentially coincide 
on the scale of the Figure (a). 

We can see that the constituent mass is almost identical to the chemical potential except 
near the canonical phase transition region. This result is the same as that of the D3/D7 
system [1] and the similar result is obtained in [25]. As mentioned in [1], this result is surprising 
because the quark density is much larger than the size of an individual quark; 



n„ 



2 3 / 2 / d \ 

-p-N f N c g eff (M q ) [ | n crit , 



(4.5) 



where g 2 eff (M q ) = \M q = g\N c M q and n 



cr it = M 3 . Structure functions [32] or a quark's 
disturbance of the adjoint fields [33] in the D3/D7 system are well studied. The mass spectrum 
of mesons in the Dp/Dq system is discussed in [35]. Motivated by their works, we assume 
that the individual quark size is ~ m gap . Combining it with the relation M/m gap ~ 0.233[12], 
we define the critical density as n cr a = M 3 at which quarks start to overlap each other. For 
the effective coupling, supergravity approximation is valid only in the region[34] 



1 « g eff « N^ 3 . 



(4.6) 



would not be negligible. 



Hence although I d/rh ) may be of order 10 



is much larger than n CT n and interactions 
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